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摘要 : 射电 观测 是 研究 太阳 活动 的 重要 探测 手段 。 我 国明 安 图 射电 频谱 日 像 仪 (MingantU 
SpEctral Radioheliograph, MUSER) 主要 用 于 研究 太阳 爆发 活动 初始 能 量 释 放 区 的 物理 过 程 ， 
其 观测 将 在 太阳 射电 成 像 开 屏 一 个 新 的 窗口 。 成 像 处 理 是 数据 处 理 的 重要 组 成 部 分 ， 如 何 提 
高 成 像 质 量 是 当前 数据 处 理 的 研究 重点 。 首 先 介 绍 了 射电 干涉 成 像 的 基本 理论 ， 随 后 分 析 了 
对 观测 得 到 的 可 见 度 数据 积分 的 必要 性 ， 细 致 讨论 了 短 时 段 可 见 度数 据 合 加 求 平均 和 长 时 段 
UV 履 盖 登 加 两 种 积分 方法 ， 并 给 出 了 完整 的 实现 。 通 过 实现 代码 与 实验 验证 ， 两 种 积分 均 
可 以 有 效 提 高 信 噪 比 ， 图 像 质量 明显 提高 。 
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射电 天 文 是 天 文学 的 一 个 分 支 ， 它 使 用 射电 望远镜 系统 在 无 线 电 波段 研究 来 自 深 空 (包括 各 类 天 
体 ) 的 射电 波 … 。 为 推动 射电 观测 技术 的 发 展 ， 填 补 国际 上 对 太阳 耀 斑 能 量 初始 释放 区 分 米 波段 高 分 
辨 射电 成 像 观 测 的 空白 和 ?1 ， 我 国 已 建成 在 0.4~15 GHz 范围 内 的 厘米 -分 米 波 日 像 仪 一 一 明 安 图 射电 
频谱 日 像 仪 。 它 可 以 在 超 沉 频带 下 同时 以 高 时 间 、 高 空间 和 高 频率 分 状 率 观测 太阳 的 动力 学 性 质 ， 探 
索 太 阳 剧 烈 活动 的 起 源 ” 。 

MUSER 天 线 阵 包含 低频 阵 ( MUSER-I1) 和 高 频 阵 (MUSER-I) 。 低 频 阵 共有 40 面 4.5 m 天 线 ， 观 
测 频率 为 400 MHz~2 GHz。 低 频 阵 的 观测 分 为 4 个 波段 ， 两 种 极 化 方式 (左旋 、 右 旋 ) ， 其 中 每 个 波 
段 包 含 16 个 通道 。 高 频 阵 共 60 TA 2 m 天 线 ， 观 测 频 率 范 围 为 2 CHz~15 GHz, Æ 33 个 波段 和 两 种 极 
化 方式 上 观测 ， 每 个 波段 包含 16 个 通道 。 高 频 阵 和 低频 阵 每 3 ms 产生 一 帧 观测 数据 (一 个 波段 、 
一 个 极 化 方向 上 16 个 通道 的 数据 ) 。 如 何 将 原始 观测 的 可 见 度数 据 经 过 一 系列 校准 和 计算 得 到 用 于 
科学 研究 的 太阳 网 像 ， 是 数据 处 理 的 关键 。 

目前 ， 国 外 用 于 太阳 射电 成 像 观测 的 设备 主要 有 日 本 野 边 山 的 太阳 射电 日 像 仪 ， 积 分 时 间 为 25 ms ° ; 
TEFL BH OK YOK BATH RAM, PAP a 1 s; 西伯 利 亚 太 阳 射 电 日 像 仪 ， 积 分 时 间 0. 336 s, X 
些 射 电 日 像 仪 积分 处 理 过 程 一 般 直 接 在 接收 机 完成 ， 因 此 时 间 分 辩 率 较 低 。 我 国 位 于 内 蒙古 明 安 图 观 
测 站 的 明 安 图 射电 频谱 日 像 仪 已 经 建成 ， 满 足 实 时 和 事后 科学 数据 处 理 的 软件 已 开发 完成 。 为 进 
一 步 提 高 成 图 质量 ， 本 文通 过 对 成 图 基本 原理 的 分 析 ， 进 一 步 提出 基于 UV 平均 积分 和 UV te PRP 
两 种 方法 提高 太阳 图 像 的 信 噪 比 。 


1 成 像 原 理 


1. 1 综合 孔径 成 像 原 理 
相关 干涉 仪 是 将 间隔 一 定 距 离 的 两 面 天 线 接收 的 天 体 辐射 分 别传 送 到 相关 需 中 进行 相关 处 理 。 由 
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干涉 成 像 原 理 可 知 ， 干 涉 仪 的 输出 啊 应 和 天 空 亮度 存在 傅 里 时 变换 关系 ， 假 设 天 空 亮度 分 布 为 了 (1， 
m), 干涉 仪 输出 为 V(u, o), ， 则 有 如 下 关系 5 ; 


Viu, %) = > > (1, mer. (1) 
将 离散 求 和 转换 成 连续 积分 可 以 得 出 如 下 公式 : 
V(u, v) = b mje S didm , (2) 


HEF, u, v 为 投影 平面 也 称 作 空间 频率 ; 1，m 为 平时 的 天 空 面 。(2) 式 表明 干涉 仪 的 输出 是 亮度 分 布 
侍 里 叶 变 换 的 一 个 (u, ov) 分量 ， 随 着 ,vw 的 不 同 可 以 得 到 不 同 空 间 频 率 分 量 的 值 ， 所 有 的 空 I 
aye ARA AY LEKA Aa E KROER EAR HF BI 8 EE, BER, ORR SE 
合 孔 径 成 像 原理 。 

对 于 任意 一 个 亮度 分 布 (1, m) ， 总 有 一 个 复 可 见 度 了 痕 数 V(u, vv) 与 之 相对 应 。 复 可 见 度 函 数 是 
连续 的 复 函 数 ， 而 天 线 阵 只 能 构成 有 限 条 基线 ， 即 只 能 获得 复 可 见 度 子 数 的 有 限 个 采样 点 。 采 样 函 数 
S(u, v)= $ ó(u -u,, v-v) ， 其 中 ,为 取 遍 所 有 (u,vw) 的 点 。 采 样 函 数 的 输出 即 为 天 线 得 到 的 真 
正 数 据 ， 对 这 些 数 据 进 行 傅 里 叶 逆 变换 得 到 脏 图 1 (1, m) U; 

I(l, m) = | s v)V(u, vye Pr”) Judy . (3) 


1.2 MUSER 积分 的 基本 原理 

当前 数据 处 理 系统 流水 线 已 完成 从 原始 观测 数据 格式 到 FITS 文件 格式 的 转换 ， 观 测 数据 校准 到 
脏 图 以 及 洁 化 的 成 图 过 程 ”。 以 高 频 阵 在 2016 年 7 月 5 日 4 时 4 分 38 秒 的 观测 数据 为 例 ， 生 成 脏 图 
(图 1(a) ) 和 洁 化 图 (图 1(b))。 


(a) (b) 


图 1 MUSER 原始 脏 图 (a) 和 洁 化 图 (b) 
Fig.1 Dirty and clean image without integration (MUSER-II, 2016-07-05 04: 04:38 (UTC) ) 


可 以 看 出 ， 对 单 帧 采样 数据 成 像 后 的 太阳 图 像 信 咯 比 很 低 ， 不 能 清楚 地 看 出 太阳 爆发 的 位 置 及 其 
轮廓 。 脏 图 洁 化 后 得 到 的 洁 图 ， 同 样 也 不 能 清楚 地 分 析出 爆发 位 置 。 为 提高 日 像 仪 重建 的 太阳 图 像 的 
信 噪 比 ， 得 到 质量 较 好 的 成 像 结 果 ， 需 要 对 采样 数据 进行 积分 处 理 ， 改 善 图 像 质量 ， 提 高 信 噪 比 ”。 
本 文 根 据 科学 研究 的 需要 ， 以 及 积分 时 段 要 求 ， 从 两 方面 分 析 ， 人 研究 了 UV 平均 积分 (可 见 度数 据 求 
和 再 求 平均 ) UV 履 善 积分 (利用 地 球 目 转 增加 UV AW) 两 种 方法 。 
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1.2.1 基于 UV 平均 的 积分 

明 安 图 射电 频谱 日 像 仪 选 择 观测 模式 为 快照 模式 ， 积 分 时 间 为 3 ms 进行 观测 。 如 此 高 的 时 间 分 
辩 率 情况 下 ， 为 了 提高 空间 分 辨 率 ， 达 到 改善 网 像 质 量 的 效 末 ， 采 取 对 3 ms 观测 得 到 的 UV 及 其 输出 
数据 平均 的 方法 实现 。 

以 高 频 阵 为 例 ， 其 中 单 帧 积分 采样 时 间 i,， 假设 对 UV 积分 时 间 为 1;,， 单 帧 采样 函数 S(u, v)， 采 
样 输出 为 V(u, v), M UV 积分 表示 如 下 : 


N 
AREF PRL : S,(u, v) = > S (us v, ) (N=—). 
i=1 


采样 输出 : k C, v) = 2 V (they v;) (Ns 


V, (u, v) 
N 
(4) RAGE BIS AE YZ YD RROGEI TEE, MR KRETE, FERE KRN n 
Jl BE KZIT EKM UV 平均 积分 。 
1.2.2 基于 UV 覆盖 增加 的 积分 
UV 窗 盖 对 成 图 起 至 关 重 要 的 作用 ， 由 于 采样 函数 S(w, >) 的 傅 里 叶 变 换 即 得 到 脏 束 ， 也 称 综合 
束 ， 而 脏 束 又 直接 影响 成 图 ， 因 此 可 通过 增加 UV BRERA eS, MN UV BA 3 种 方式 ， 
(1) 增 加 阵列 望远镜 个 数 ;，(2) 增 加 子 孔径 有 效 半径 ;(3) 利 用 地 球 自转 "” 。 增 加 阵列 望远镜 个 数 对 
成 本 有 直接 需求 ， 在 成 本 固定 的 情况 下 不 考虑 ， 增 加 子 孔径 有 效 半 径 可 能 导致 采样 率 不 足 。 在 现 有 成 
本 下 通过 地 球 自转 增加 UV 履 盖 是 一 种 非常 直接 可 行 的 方法 。 
随 着 地 球 自转 UV 发 生 改 变 ， 而 短 时 间 内 UV 的 变化 很 小 ， 因 此 假设 每 间隔 i 增加 UV 覆盖 ， 
总 积分 时 间 t, ， 则 增加 UV 覆盖 后 的 采样 函数 和 采样 输出 表示 为 


— ` x: y t, 
ARIE PRB: S, (u, v) = DSi v) L ) . 
j=l 


S (u, v . 
UV 平均 脏 图 : pU, = JJ — a — (w=). (4) 


采样 输出 : V u, o) = Vw 6) (N=). 


UV 覆盖 增加 脏 图 : O, m) = fS, Cu, v) V, (u, ve mem) dud , (5) 


t, 


由 (5) 式 可 知 ， 在 间隔 ww 增加 UV Bimi, ARISA] 可 以 增加 - 倍 的 UV 数据 点 ,假设 有 
interval 
n 面 天 线 构成 的 天 线 阵 列 ， 可 以 得 到 n(n-1)/2 个 数据 采样 点 ， 通 过 增加 UV bith, 采样 数据 点 总 共 
t 1 ‘ ep as ` EL: ¿H -L- 
al (EOD) S, EAER REOR 


[inal 
上 述 两 种 方式 均 可 以 提高 成 图 的 质量 ， 将 两 种 方法 混合 到 一 起 同样 可 以 提升 图 像 质量 ， 由 此 提出 
先进 行 UV 平均 积分 再 增加 UV 履 盖 积分 的 方法 。 
1.2.3. 基于 UV 平均 和 增加 覆盖 的 混合 积分 

假设 UV 积分 时 间 为 1;,， 可 以 知道 在 该 段 时 间 内 的 采样 函数 为 S Cus 0), 采样 输出 V, (ui, s), 
每 间隔 上 ,增加 UV 覆盖 ， 总 积分 时 间 ， ， 混 合 积分 采样 函数 及 其 采样 输出 可 以 表示 为 


M 
vy È Sw, vy) 


— P. sp: t; t, 
RAE PRAT: S a (u, v) = > (M=—, N= fa 
j=l M t, Lint al 
M 
v 之 Val Wy, vy) ! ) 
采样 输出 : l sm, v) T > - (M=—, N= - ) 。 
j=l M t bint al 
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HEAL: I, m) = [[ Suu, v)V,,(u, veo” dudo . (6) 


t, 


在 总 时 间 t, 内 UV 覆盖 增加 了 音 ， 每 组 增加 的 UV 禾 新 是 在 时 间 段 tj 对 UV Fea Feeney. 


interval 
2 ”积分 过 程 实现 


明 安 图 射电 频谱 日 像 仪 的 成 图 过 程 : 目标 源 放射 出 天 体 辐射 ， 然 后 干涉 仪 对 其 进行 采样 ， 采 样 得 
到 的 原始 数据 经 过 预 处 理 相 关 校 准 … 、 对 应 极 化 和 通道 ， 之 后 进行 积分 处 理 ， 再 网 格 化 传 里 叶 逆 变 
换 得 到 脏 图 ， 最 后 经 过 洁 化 算法 处 理 得 到 洁 图 ， 过 程 如 图 2. 


采样 数据 积分 处 理 


Clean、 Steer 


图 2 MUSER 积分 成 图 过 程 
Fig.2 The procedure of MUSER integration 
2.1 积分 操作 流程 
明 安 图 射电 频谱 日 像 仪 的 每 3 ms 接收 一 帧 观测 数据 ， 对 这 些 数据 进行 积分 操作 得 到 积分 数据 ， 
积分 数据 是 观测 得 到 并 且 经 过 预 处 理 相关 校准 、 对 应 极 化 和 通道 的 可 见 度数 据 。UV 平均 积分 具体 操 
作 流 程 如 图 3。 


完整 数据 帧 完整 数据 帧 完整 数据 帧 平均 后 同样 
完整 数据 帧 完整 数据 帧 Zü J i Ne Md 


2.0~2.4 GHz 


14.8~15.2 GHz 14.8~15.2 GHz 


, — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — 


图 3 可 见 度数 据 平均 的 计算 


Fig.3 Calculation of the average of visibilities 


高 频 阵 中 每 一 帧 完整 的 数据 包含 33 小 帧 。UV 积分 对 连续 的 每 一 帧 数据 中 UU, VV 以 及 采样 得 
到 的 实 部 、 虚 部 数据 进行 累加 再 求 平 均 ， 得 到 的 UV 平均 积分 数据 量 和 单 帧 大 小 数量 一 样 多 ， 例 如 在 
高 频 阵 中 采样 函数 UV 点 为 3 540 个 采样 得 到 的 实 部 虚 部 3 540 个 ， 则 经 过 UV 平均 积分 得 到 的 数据 同 
样 是 3 540 个 UV 点 及 其 采样 值 。 


82 


图 像 的 质量 ， 提 高 信 噪 比 。 


采样 数据 进行 平均 ， 之 后 每 隔 10 min 增加 UV 22 m, 
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实现 UV 履 盖 积分 具体 操作 流程 如 图 4。 对 高 频 阵 中 数据 进行 UV 履 盖 积分 同样 用 完整 的 数据 帧 
积分 ,不 同 的 是 由 于 短 时 间 内 采样 函数 及 其 采样 值 都 变化 很 小 ， 如 果 用 连续 的 帧 增加 UV 48 w Y 2828 
有 末 不 好 ， 为 此 需要 隔 一 段 时 间 增 加 UV Zeb ait. UV 黎 盖 积分 得 到 的 是 一 段 时 间 间 隔 完 整数 据 帧 三 加 的 
数据 ， 同 样 对 数据 UU, VV 及 其 采样 值 操 作 得 到 的 数据 量 成 倍增 加 ， 这 种 操作 能 极 大 改善 稀 踊 采样 


完整 数据 帧 


overlap 


14.4~14.8 GHz 
14.8~15.2 GHz 


完整 数据 帧 


3.2~3.0 GHz 


14.4~14.8 GHz 


14.8~15.2 GHz 


overlap...overlap 


2.0~2.4 GHz | 
2.4~2.8 GHz | 
2.8~3.2 GHz 


14.4~14.8 GHz 
14.8~15.2 GHz 


图 4 UV RRR RAE 
Fig.4 Procedure of UV coverage overlap 


U8 tr ee JII 
数据 帧 


2.0~2.4 GHz 


14.8~15.2 GHz 


以 积分 60 min 的 数据 为 例 ， 积 分 后 增加 UV Hiatt, SRC pR APO H, XH 02 RE Te au Es] E Dr t 
及 其 信 噪 比 起 至 关 重 要 的 作用 。 高 频 阵 单 帧 的 UV # w 4 S BOE 1770 4, SESE SRA 3 540 个 数 
据 ， 得 到 一 张 2 560 x 2560 像素 的 图 像 ， 窗 盖 率 仅仅 0. 054% 。 而 对 数据 总 共 积 分 60 min 间隔 10 min 
AD UV A maA 10620 个 数据 点 ， 共 斩 之 后 数据 点 21240 4", AKESA 0.324%. HR 60 ms 
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图 5 3ms 单 帧 UV Re 
Fig.5 UV Coverage for a single frame (3mas ) 
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得 到 的 UV 覆盖 图 像 如 图 $、 图 6。 
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图 6 60min 积分 UV BA 


Fig. 6 


UV Coverage for 1 hour’s integration 
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2.2 积分 代码 实现 


为 了 实现 UV 平均 积分 、UV ERIAK UV 平均 和 UV 和 窗 盖 混合 积分 ， 采用 Python 编写 ， 同 
时 为 了 简化 程序 ， 将 3 种 模式 的 积分 写 人 一 个 接口 ， 提 供 模式 的 选择 ， 程 序 接口 形式 如 下 : 


def clean_integration_R (self, sub_ARRAY, is_loop_mode, start_time, end_time, TASK TYPE, time_average, 
time_interval, BAND, CHANNEL, PLOT ME, WRITE_FITS, P_ANGLE, DEBUG, outdir) : 


主要 参数 : sub_ARRAY 为 MUSER 期 数 的 选择 MUSER-I sk MUSER-II, start_time 和 end_time 为 
开始 积分 时 间 和 结束 积分 时 间 ，TASK_TYPE 为 积分 模式 的 选择 包含 UV 平均 积分 、UV 覆盖 积分 、 
UV 平均 和 UV IEA AR, time_average 为 UV 平均 积分 的 积分 时 间 ，time_interval 为 UV 72 te FA) 
时 间 间 隔 即 间隔 多 长 时 间 闪 加 一 次 UV 和 覆盖。 在 系统 中 的 运行 情况 如 图 7。 代 码 编写 可 以 适用 于 数据 
采集 的 任何 时 间 段 积分 模式 。 


Muser<1>: inp integrationuvfits 
一 一 一 一 一 一 一 一 > inp(integrationuvfits) 
# integrationuvfits :: Integration FITS file for MUSER 


subarray = 1 # MUSER I or II (1, 2) 

is_loop_mode = True # LOOP Mode or not(True, False) 

start_time = ie # start time (start_time = "2016-07-05 
# 12:04:38 463") 

end_time - # end time (end_time = "2016-07-05 
# 12:30:00") 

task_type = sl # task type: average, interval or 
# mixture 

time_average = .. # time span of average computing 

time_interval = None # time_interval 

debug = 0 # Display increasingly verbose debug 
# messages (0, 1) 


图 7 MUSER 系统 运行 积分 命令 


Fig.7 Integration task in MUSER data processing System 


积分 实例 以 混合 积分 模式 为 例 ， 进 行 10 min, 20 min, 30 min, 40 min, 50 min 以 及 60 min 的 数据 
做 积分 处 理 ， 对 60 ms 的 数据 进行 UV 平均 积分 ， 之 后 每 隔 10 min 增加 UV 覆盖 。 脏 图 结果 如 图 8。 

图 8 从 左 到 右 从 上 到 下 分 别 为 积分 10，20，30，…，60 min 之 后 的 脏 图 ， 可 以 很 明显 地 看 出 从 第 
1 张 到 最 后 一 张 太阳 轮廓 略 清晰 。 随 着 积分 时 间 的 增加 得 到 的 脏 图 质量 越 来 越 好 ， 图 像 越 来 越 清 晰 ， 
同时 可 以 清楚 地 看 到 太阳 的 轮廓 (图 中 间 类 似 圆 形 ) 以 及 太阳 爆发 的 亮点 位 置 ， 很 明显 积分 达到 提升 
图 像 质量 的 效果 。 经 过 洁 化 处 理 得 到 的 洁 图 如 图 9。 

图 9 从 左 到 右 从 上 到 下 分 别 为 积分 10，20，30，…，60 min 之 后 的 洁 图 ， 可 以 很 明显 地 看 出 ， 随 
着 积分 时 间 的 增加 ， 得 到 的 洁 化 图 像 越 来 越 清 晰 ， 同 时 可 以 清楚 地 看 到 太阳 的 轮廓 以 及 太阳 爆发 的 位 
置 亮点 ， 从 一 开始 多 处 亮点 爆发 到 最 后 一 张 只 有 一 处 非常 明显 的 爆发 点 ， 信 品 比 有 很 大 提升 。 

图 10 为 2016 年 7 月 5 日 04:04:38UT， 射 电 日 像 仪 高 频 阵 在 4. 187 8 GHz 上 用 60 min 积分 时 间 
观测 得 到 的 太阳 射电 图 像 结果 。 与 日 本 野 边 山 射电 日 像 仪 2016 年 7 月 5 日 04:00:02UT 在 17 GHz 上 
观测 得 到 的 太阳 图 像 比较 ， 高 频 阵 得 到 的 太阳 像 与 野 边 山 射 电 日 像 仪 观测 得 到 的 太阳 像 结果 基 本 一 
致 ， 可 见 积分 方法 的 实现 结果 比较 可 靠 。 
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Fig.8 Dirty images of mixed integration 
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图 9 混合 积分 洁 图 


Fig.9 Clean images of mixed integration 
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图 10 (a) MUSER-II 60 min 2265 AMER; (b) 野 边 山 射 电 日 像 仪 观测 的 太阳 像 
Fig. 10 (a) One huor's integration solar images of MUSER-II; 


(b) Solar observation images of the Nobeyama Radioheliograph 
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通过 对 采样 数据 积分 处 理 ， 以 每 60 ms UV 积分 和 间隔 10 min 增加 UV 覆盖 混合 积分 方法 为 例 ， 
可 以 看 出 成 图 效果 非常 明显 ， 经 过 一 系列 数据 人 处理 计算 得 到 的 脏 图 可 以 清楚 地 看 出 太阳 边缘 轮廓 及 其 
爆发 位 置 ， 信 噪 比 有 了 很 大 提升 。 并 且 随 着 积分 时 间 的 增加 ， 得 到 的 图 像 质量 也 有 所 提升 。 这 一 积分 
思路 对 于 获得 质量 较 优 的 MUSER 观测 结果 有 较 好 的 作用 。 目 前 图 像 处 理 结果 得 到 很 大 提升 ， 但 在 成 
图 前 需要 读 取 UVFITS 数据 文件 ， 当 前 的 UVFITS 文件 生成 时 间 较 慢 ， 急 需 加 快 生成 的 效率 。 后 续 积 
分 工作 可 以 在 积分 时 间 上 进行 探索 ， 利 用 评价 指标 寻找 如 何 积分 才能 达到 相对 比较 好 的 效果 ， 给 出 
UV 积分 时 长 、UV 履 盖 积分 时 间 间 隔 以 及 总 共 60 min 积分 时 间 等 参数 的 量化 值 。 
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Integral Method and Implementation of MUSER Visibility Data 
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Abstract; Radio observation is the most important detection means to research the solar activity. The 
MingantU SpEctral Radioheliograph (MUSER) is mainly used to research the physical process in the initial 
energy-releasing area during the solar activity and its observation will open a new window in solar radio 
imaging. Imaging processing is an important part of MUSER data processing and improving the image quality is 
the most important research of data processing. Firstly, we introduce the basic theory of radio interference 
imaging. Secondly, we analyze the necessity of visibility data integration obtained from the MUSER 
observations. Then, we discuss in detail the two kinds of integration method which are short-term visibility data 
superimposed averaging and long-term UV overlay stacking and give out a complete realization. Through the 
method of verification and the validity of completion, the two kinds of integration method can effectively 
improve the signal-to-noise ratio and the quality of images. 


Key words: MUSER; UV integration; UV coverage integration; Aperture synthesis imaging 


